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Abstract 

CO ' 

| The theory relevant to the study of matter in equilibrium with 

the radiation field is thermal quantum electrodynamics (TQED). We 
present a formulation of the theory, suitable for non relativistic fluids, 
\ based on a joint functional integral representation of matter and field 

ON ' variables. In this formalism cluster expansion techniques of classical 

statistical mechanics become operative. They provide an alternative 
| to the usual Feynman diagrammatics in many-body problems which is 

not perturbative with respect to the coupling constant. As an applica- 
tion we show that the effective Coulomb interaction between quantum 
charges is partially screened by thermalized photons at large distances. 
More precisely one observes an exact cancellation of the dipolar elec- 
tric part of the interaction, so that the asymptotic particle density 
correlation is now determined by relativistic effects. It has still the 
q-( r D decay typical for quantum charges, but with an amplitude strongly 

reduced by a relativistic factor. 

•i-H . 

^ ! 1 Introduction 

A precise and complete description of equilibrium states of non rel- 
ativistic quantum charges interacting via the static Coulomb poten- 
tial has been thoroughly developed in recent years in the low density 
regime [l]-[5]. This description relies on the use of the Feynman-Kac 
path integral representation of the thermal Gibbs weight allowing for 
a classical-like analysis of thermodynamic potentials and particle cor- 
relations. Essentially, quantum point charges are mapped onto a set of 
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closed Brownian paths (loops) whose random shapes account for the 
quantum fluctuations. Techniques of classical statistical mechanics be- 
come available in the auxiliary phase space of loops, in particular the 
method of cluster expansion (Mayer graphs) . The latter is particularly 
suited to calculations in dilute systems, where the small parameter is 
the density. 

Low density expansions of the pressure are performed up to the or- 
der p 5 / 3 [1], exact asymptotics of particle correlations are determined 
in [2]. Phases with atomic or molecular recombination can also be 
conveniently studied, e.g. the equation of state [3] and the van der 
Waals forces [4] in the Saha regime, as well as the dielectric response 
of an atomic gas [5] (see [6] , [7] for reviews and additional references) . 
However, none of these works take into account the coupling of the 
charges to the radiation field which is responsible for both effective 
magnetic interactions (Lorentz forces) and retardation effects. The 
purpose of this paper is to show how the above formalism and tech- 
niques can be generalized when matter is thermalized with the quan- 
tized electromagnetic field. It is an extension of [8] (hereafter referred 
to as I) where the field was considered as classical. When the field is 
quantized in the transverse gauge, it is appropriate to represent the 
Gibbs weight by means of the bosonic functional integral based on the 
coherent state representation of photon states. In this way the quan- 
tum field is mapped onto a set of classical-like random electromagnetic 
fields with (imaginary) time dependent amplitudes. Since the energy 
of the free field is quadratic in the field amplitudes, the latter are dis- 
tributed with Gaussian statistics. At this stage, quantum charges can, 
as in [l]-[5], be put into correspondence with Brownian charged loops 
with the aid of the Feynman-Kac-Ito formula. The coupling to the 
field appears as the flux of the magnetic field accross the loops. Thus 
TQED becomes isomorphic to a system of random charged wires (the 
loops) experiencing a random magnetic field. The calculation rules 
are entirely defined by the covariances of the processes associated to 
the loops and to the field amplitudes, together with the use of Wick's 
theorem. In this setting, the cluster (Mayer or virial) expansions of 
classical statistical mechanics can again be put at work providing an 
alternative to the standard TQED Feynman graph calculations, which 
is not perturbative with respect to the coupling constant (namely, the 
electric charge). The method is particularly adapted to study equi- 
librium phases of plasmas and recombination processes in presence of 
the electromagnetic field at moderate density . 
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In Section 2 we describe the actual system consisting of non rel- 
ativistic charges interacting with the photon field. In order to make 
sense mathematically and physically, the model requires a high en- 
ergy cut-off defined by ^k cut = ™c 2 to eliminate photons that are 
more energetic than the rest mass energy of a particle of typical mass 
fh (H is the Planck constant, c the speed of light and uj^ = ck the 
photon frequency for the wave number k, k = |k|). This gives a typi- 
cal wave number cut-off k cn t = fhc/h with corresponding wave length 
A C ut = h/fhc (see e.g. [9], Chap. 3, for a discussion of this point). 
High energy processes, such as pair creation or annihilation, demand 
for the use of the relativistic wave equation (Klein-Gordon or Dirac). 
They are not taken into account in this model whose predictions only 
make therefore sense for distances r 3> A cut . 

The construction of the relevant functional representations are re- 
called in Section 3 for the field and in Section 4 for the particles. Since 
the subject is well developed elsewhere we merely present the main 
structure in a perspective adapted to our purposes (see references in 
Section 3). The thermalized photon field involves the typical energy 
^k ph = ^cfcph = ksT = (5~ l , with corresponding wave length A p h = 
Phc, called the thermal length of the photon (T is the temperature and 
fee the Boltzmann constant). On the other hand, the mean kinetic en- 
ergy of a nonrelativistic particle £k mat = (^mat) 2 /2fTi = /c#T defines 
the de Broglie thermal wave length of the particle A mat = h^/ (3/fh. 
To be consistent with non relativistic particle motion we must im- 
pose that the thermal energy imparted to the particle in the form 
of kinetic energy is much lower than its rest mass energy, namely 
£ k mat =k B T < rhc 2 , implying 

A cut = < A mat <C Aph = V /3mc 2 A mat (1) 

y pmc 2 

where fifhc 2 > 1 is a dimensionless relativistic parameter. There- 
fore, when the field is quantized, we have to distinguish two different 
regimes at large distance r 

Amat < Ap h < r (2) 

or 

Amat < r < Aph . (3) 

In addition to the quantum lengths, there are typical classical lengths 
such as the interparticle distance a = p -1 / 3 (p the density) and the 
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Debye screening length Ad- The latter do not enter explicitly in our 
subsequent analysis because the regime (2) of main interest in this 
paper deals with distances r far beyond A mat , a and Xp. For instance, 
in an electrolyte the lengths A mat , a and Xd are of the same order of 
magnitude (~ 10 -10 — 10~ 9 m) but they are all much smaller than 
A p h ~ 10~ 5 m (see concluding remarks). We shall only require that 
the density is low enough for the system to be in a fluid phase so that 
we can apply the standard methods of statistical mechanics (cluster 
expansions). 

In Section 5, we determine the effective potential between loops 
arising when the field degrees of freedom have been integrated out. 
This can easily be done by a Gaussian integration, as in paper I. In- 
deed, a simple structure shows up from the fact that in the functional 
integral representation the coupling of matter to the field amplitudes 
occurs linearly in a phase factor (in contrast to the original quantum 
Hamiltonian which has a coupling quadratic in the creation and anni- 
hilation operators). Then the whole effect of the field is contained in 
an effective potential depending on A ma t and A p h that can be viewed 
as a current-current interaction between pairs of loops (Formula (66) 
in Section 5). 

We use these results in Section 6 to find the behaviour of the 
particle correlations in both regimes (2) and (3). Equipped with the 
Coulomb potential and this new effective field-induced potential, all 
standard rules of classical statistical mechanics can be applied to the 
calculation of particle correlations (some care has to be exercised with 
the computation rules for stochastic integrals, see appendix A). It is 
seen that the large distance behaviour of the correlation is determined 
by the square of dipoles fluctuations, the total dipole of a loop having 
a part due to its charge and a part due to its current. This leads to 
a generic r -6 decay of the correlation. Now a striking phenomenon 
occurs in case (2) above: namely the screening of the dominant part 
of the Coulomb interaction by thermalized photons. When r 3> A p h, 
the transverse field has a contribution that exactly cancels the dipolar 
electric part of the loop fluctuations. Only current fluctuations of the 
loops are left, which cannot be screened. In this regime, the correlation 
still has a r -6 decay, but with a relativistic prefactor {[5mc 2 )~ 2 . 

In paper I, we have argued that large distances are controlled by 
small wave numbers of the radiation field and the latter can therefore 
be treated classically. This apparently sensible argument proves to be 
incorrect in the sense that it does not predict the aforesaid Coulombic 
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cancellation which results of a subtle conspiracy between the Planck 
constants of field and matter. It might be inconsistent, in the trans- 
verse gauge, to make a classical approximation for the radiation part of 
the field only. Approximations should be made in a fully gauge invari- 
ant manner. Note however that, once the cancellation has been taken 
into account, the theory of paper I correctly predicts the remaining 
correlation tail induced by the current fluctuations. 

In the regime (3), the radiation field has essentially no incidence 
on the decay of the particle correlations and one recovers the purely 
Coulombic tail due to electrical dipole flucuations as the dominant 
contribution, plus terms vanishing as r/A p h — ► 0. More generally, all 
results of [l]-[5] are expected to remain valid in this regime up to tiny 
relativistic corrections. 

Other applications for which the present formalism will be relevant 
are suggested in the concluding remarks. 

2 The model 

The non relativistic QED model consists of non relativistic quantum 
charges (electrons, nuclei, ions) with masses m 7 and charges e 7 . They 
obey the appropriate Bose or Fermi statistics and interact with the 
quantum electromagnetic field, the latter being relativistic by nature. 
The index 7 labels the S different species and runs from 1 to S. The 
particles are confined in a box A £ R 3 of linear size L whereas the 
field itself is enclosed in a large box K with sides of length R, R^> L. 
The Hamiltonian of the total finite volume system reads in Gaussian 
units 

H l ,r - }^ ^ ^ lr- -r-l ^ ^t(7i, Ti > + H o 

i=l 7l i<j 1 * - 71 i=l 

(4) 

The sums run on all particles with position r-j, momentum pi and 
species index 7$, i=l,. . . ,N , V c ^t{^i^i) comprises a possible external 
potential plus a steep wall potential that confines the particles in A. 
The latter can eventually be taken infinitely steep at the wall's position 
implying Dirichlet boundary conditions on the particle wave functions 
at the boundaries of A. 

The electromagnetic field is written in the Coulomb (or transverse) 
gauge so that the vector potential A(r) is divergence free and H™ d 
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is the Hamiltonian of the free radiation field. We impose periodic 
boundary conditions on the faces of the large box K. Expanding 
A(r) and the free photon energy in the plane wave modes k = 

( ^2^2™, )gives 

Mr) = (^) V2 E^W^(4Ae" ikr + a k Ae^) (5) 

kA 

H^ d = ^2hw k al x a kX , (6) 

kA 

where aj. , a k x are the creation and annihilation operators for photons 
in the mode (kA) with commutation relations [a^x, a^y] = ^AA'^kk'i 
©kA (A = 1,2) are two unit polarization vectors orthogonal to k and 
uo k = ck, k = |k|. In (5) g(k) is a real spherically symmetric smooth 
form factor that takes care of the ultraviolet divergencies. It obeys 
g(0) = 1 and is supposed to decay rapidly beyond the characteristic 
wave number fc cut = fhc/h. Note that in (4) we have included neither 
the Pauli coupling — /x • B(r) of the electronic spin with the magnetic 
field B(r) = V A A(r) (ji = {eh/ Am e c)a is the magnetic moment 
of the electron, <x are the Pauli matrices) nor the nuclear hyperfine 
interaction (see comments in the concluding remarks). It is known 
that the Hamiltonian (4) is ^/-stable [10] for a finite ultraviolet cutoff 
k^ 1 < oo, namely Hl,,r possesses an extensive lower bound propor- 
tional to the total number of particles (for a review of ii-stability in 
non relativistic QED, see [11]). 

We are interested in the situation in which matter and photons are 
in thermal equilibrium at the same temperature T. The total partition 
function associated with (4) 

Z LyR = Tr e-^.« (7) 

is obtained by carrying out the trace Tr = Tr ma tTr ra d of the total 
Gibbs weight over particles' and the field's degrees of freedom, namely 
over the particle wave functions with appropriate quantum statistics 
and the Fock states of the photons. The corresponding free energy 
density in the thermodynamic limit will be defined by extending to 
infinity first the field region K and then the box |A| containing the 
charges. Thus the excess free energy relative to that of the free radi- 
ation field is 

/ = -k B T lim -L lim (In Z L)R - In Z$) , (8) 
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where = Tr ra( j exp (— (5Hq^) is the partition function of the free 
field. A lower bound for / has been established in [12], but at the 
moment, to our knowledge, a complete proof of the existence of the 
thermodynamic limit has not yet been provided. Nevertheless we shall 
assume that the quantities of interest in this paper have a well-behaved 
thermodynamic limit. 

As in I, we shall be concerned in the sequel with the partial average 

Tr jp"^i 1 « 
[e-^-Wt = Trrad ! rad 0) 

giving the (non normalized) statistical distribution of matter obtained 
by averaging on the degrees of freedom of the radiation field. The 
corresponding normalized reduced density matrix is 4 

Tr rad e-^* [e-^^Ut nn , 

PL ' R ~ z L , R ~ M^U ' ( 0) 

It will be convenient to single out in B.l,r the free radiation part 
writing 

H R , L = H a + H™*, (11) 
N f D . _ ^LA(r )) 2 

i=i 



where 



N e e N 



j y -\ + Z^VextilhTi) (13) 
i<j Tj »=1 

is the total potential energy. 



3 Functional integral representation of 
the field 

If the field is treated classically (namely the creation and annihilation 
operators are replaced by c-number amplitudes) it is immediately seen 

4 Here the notation is slightly different from paper I where Pl,r (1.5) designates the 
field averaged quantity (9). 
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that the free field distribution factorizes in the total Gibbs weight as 
exp (— PHr : l) = exp (— /3i?Q ad ) exp (— [5Hj±). Thus the partial trace 
(9) reduces to integrals with a Gaussian weight since the free radiation 
part exp (— /3i?Q ad ) is Gaussian in the field amplitudes, a fact that was 
exploited in I. 

If the field is quantized, it is first of all necessary to represent 
the electromagnetic field by c-functions in the total Gibbs weight 
exp (—PHr : l)- This can be achieved by means of the standard func- 
tional integral for bosonic quantum field [13], [14]. We briefly recall 
its construction. First, one considers the coherent states associated to 
the field modes 



oo 



( t V 

Ka> = £ ^ x y> |0) = e^A| >, «kA|«kA) = a k A|a kA ) 

771. 

(14) 

They have scalar products 

(ctkAK'A') = e Q kA<'A' (15) 
and the closure relation reads 

>kA|2 KA><a k A| = 1 • (16) 



/ 



d 2 a k A 

e |« k A/\«kA| = J- • 

d 2 C 



7T 



We denote a = {a k A}kA, \a) = Uk\ \ a k\), da = ElkA — aa = 
SkA a kA«kA' e ^ c -' an< i introduce the infinite product representation 
e-f> H = lim M ^oo (l - ^^) M where H = H L , R = H A + ^ ad is 
the total Hamiltonian operator (11). Using this representation and 
inserting M — 1 closure relations one can write the following coherent 
state matrix element as 



(a\e~ pH \a) = lim 



M^oo 



M-l r 

J] date-*!* 



(17) 



x (a\ (l - &H) |a M -i) • • • (l - &h) |a*-i> ■ ■ ■ <<*i| (l - &H) \a) 

As a first step we consider the partial coherent state matrix element 
{(Xi\e~P H \aii-i), which is still an operator acting on the Hilbert space of 
the particle states. Its evaluation is achieved by putting H in normal 
order. Using (15), this yields 



(<*l\ (] 



&H) |a,_!) = e a >-i(i _ JjLH(a*,oti-\)), (18) 
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where H{a1 , tty_i) depends on the complex amplitudes a according to 
the normal order form of H. Prom (11), (12), one finds 



fl"(a,*,a,_i) =H A (alai- 1 )+D N + H I ad (alai- 1 ) 



N 



# A «,a*-i) = ^ 



i=i 



(Pi ~ ^fMr^o^a^Y 

2m 7i 



+ ^ P ot(ri,7i, ■ ■ ■ ,*N, 
(19) 



where the vector potential A^^a^an) has the same form as in 
(5) with the operators aj^,akA replaced by the complex amplitudes 
a z*kA' a /-i,kA, and likewise for fl5 ad (a i *,ai_i). The constant 



D 



A' 



c m n 



?3 



R 3 



g 2 (k) 
fc 



(20) 



arises when putting (A(r-j)) in normal order. Inserting (18) in (17) 
yields 



{a\e-P H \a) 



lim 

M^oo 



M-l 

i=i 



x ( 1 - ^H{a*,a M -i. 



^fr(a,*,a,_i; 



a; 



(21) 
fl-(aj,a)) 



One introduces the formal functional integral as usual by interpreting 
a J,kA = a kA(]^) as the value at r = of a closed trajectory akA(7~) 
in the complex plane, akA(O) = ctkA(l) = «kA- The parameter r, 
< r < 1, is a dimensionless imaginary time. In the limit M — > oo 
the product of infinitesimal evolutions in (21) tends to the time ordered 
propagator 

e -f3D NT [ e -f3f ( l drH(a*(r+r,),a(r))j = 

Q -PD N e -f3ti dr^( a «(r+ f) ), a (r)) r j e ^/ 1 drH A ( a «(r+,),a(r)) j ^ ^ _^ Q + 

(22) 

The imaginary time ordering T is necessary because although the field 
amplitudes ot(r) are now c-functions, the H(a*(r + 77), a (r)) are still 
operators acting on the space of particle wave functions and there- 
fore they do not commute for different times. However, the free field 
part #Q ad (a:*(T + 77), oc(t)) commutes with the matter dependent part 
Ha{<x*(t + v)i a i T )) (19) an d can b e factorized out of the T-product 
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according to the second line of (22). The rj — > + prescription means 
that, as a result of the normal order, the amplitudes correponding to 
the creation operators a*(T + 7/) have to be evaluated in (22) at times 
infinitesimaly larger than those corresponding to the annihilation op- 
erators a(r) (see (21)). Finally, (21) can be written in the condensed 
form of a path integral 



(a\e~ pH \a) = c 



-f3D„ 



lim 



o(l)=o 



d[tt(-)] e" J <fr(«^)£a(T)+/9tf5*(«(T))) 



a(0)=a 

r[ e ^/o dTff A(«W)] 



(23) 



where the bracket [■ ■ ■ ]^ indicates that the amplitudes a* in (23) have 
to be evaluated at the time r + rj. The partial Gibbs distribution (9) 
is obtained by integrating the matrix element (23) on dot and then 
dividing it by the partition function of the free field 



[e-^-U = ^/cfa e -W a <a|e-^|a) . 



More generally, the factor 

- £ dr (a* (r) £a(r)+f1H™* (o(t))) 



(24) 



(25) 



in (23) provides a Gaussian (free) weight on the space of time-dependent 
complex field amplitudes o>(t). If F(a(-)) is a functional of these am- 
plitudes, we will denote its average with respect to the distribution 
(25) by 



< F(Of(-)) >rad= 



j lim 



(26) 



Here the integral runs over all possible closed paths, setting J Da ■ 



It is well known that this Gaussian inte- 



gral is characterized by the covariance [13] 

(a k A (t)o&/ A' ( r '))rad = ^AA'^kk'C (fc, T - t') 
(a k A(r)Q k /A'(r / ))rad = (akA( r )«k'A'( T/ ))rad = 



(27) 
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with 

C{k,r- t') = e - f5nw ^ T - T '\6{T - r')(n k + 1) + 6{t' - r)n k ], r ± r' 

(28) 

C(k,0)=n k , t = t' (29) 

and 

n k = (e^ k - I)" 1 (30) 

is the Planck distribution (6 is the Haevyside step function). The 
function C(k,T — T r ) is discontinuous at r = r' with the value C(&, 0) = 
n k as a consequence of the normal order prescription rj — > 0+. 

Functional integrals (26) of the paths ck(t) are in principle entirely 
determined by application of Wick's theorem and use of the covariance 
(27). In particular, using the representation (26), the effective partial 
thermal weight (9) of matter when the field degrees of freedom have 
been traced out can now be written as 

[e -W L , R]mat = e -PD N <T[e -/J/ 1 drH A («(r))] ^ . (31) 

This will be the starting point of our investigation of the particle 
correlations in presence of the field in Section 5. 

4 Functional representation of the par- 
ticles 

We now come to the functional integral representation of the matter 
degrees of freedom. One notes that the operator T[e - ^o dr ^A(<*(i-))j 
in (31) is the propagator on the space of particle wave functions asso- 
ciated to the time dependent Hamiltonian Ha(oi(t)) where the vector 
potential has been replaced by its non-operatorial classical form 

A(r,a(r))= (^y /2 ^ 5 (fc)-^=(a kA (T)e- k ' r + a kA (r)e^) . 

(32) 

The time dependence is introduced by the amplitudes a k A(r), which 
are random functions distributed by the Gaussian weight (25) of the 
bosonic functional integral. However, for a fixed function 
ot(r), < r < 1, Ha(oi(t)) can be viewed as the Hamiltonian of the 
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particle system submitted to the external vector potential (32). In this 
situation one can apply the Feynman-Kac-Ito formula [16] to represent 
the configurational matrix element of T[e - ^o dr ^A (<*(») j _ 

For a single particle of mass m and charge e in a scalar potential 
V ext {r) and time dependent vector potential A(r, s), we first recall 
that this matrix element reads [15], [16], [17] 

v2 i\ , . x 3 /2 



(r|T exp \ -/3 / dr 



(p-§A(r,r))' 



2m 



+ V 



cxt / 



|r> = 



2vrA 2 J 



D(|) 



exp 



-P 



[dr V cxt (r + \i(r)) 
Jo 



/W(r)-A(r + ^(r),r) 
Jo 

(33) 



Here£(r), < r < 1, £(0) = £(1) = 0, is a closed dimensionless Brow- 
nian path and D(£) is the corresponding conditional Wiener measure 
normalized to 1. This measure is Gaussian, formally written as 

2 . 



D(0 = exp 



It has zero mean and covariance 





d£(r) 




dr 



)dK(-)] • 



|D(0 ^(r)r(r) = V(mm(r, r') - tt>) , 



(34) 



(35) 



where £ M (r) are the Cartesian coordinates of ^(r). In this represen- 
tation a quantum point charge looks like a classical charged closed 
filament T = (r, £) located at r and with a random shape £(r), < 
r < 1, the latter having a spatial extension given by the thermal de 
Broglie length A = fry 7 j3/m (the quantum fluctuation) . The magnetic 
phase in (33) is a stochastic line integral: it is the flux of the mag- 
netic field across the closed filament. The correct interpretation of this 
stochastic integral is given by the rule of the middle point, namely, 
the integral on a small element of line x — x' is defined by 

d£-f(0 = (x-x')-f f^^y x-x'^0. (36) 



We shall stick to this rule when performing explicit calculations 5 . If 
there is no field, the generalisation of the Feynman-Kac formula to the 



5 We find it convenient to apply the middle point rule because it correctly represents 
the quantum mechanical Gibbs weight in presence of a vector potential (divergenceless 
or not) [16]. Although we shall not use the Ito prescription we keep the terminology of 
Feynman-Kac-Ito formula. 
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many particle system including quantum statistics has been presented 
in a number of works, see e.g. [2], [7], [18]. When the field is present, 
the analysis presented in the above works can be reproduced without 
changes, the only difference being the inclusion of the additional phase 
factor corresponding to the vector potential (see [19] in the case of 
a uniform magnetic field). We give here merely the basic formulae 
resulting from these generalisations. 

Filaments T = (r, £(r), < r < 1) associated to single quantum 
particles are generalized to Brownian loops 

£=(r, 7 ,g,X(r)), < r < q . (37) 

The g-loop C consists again in a closed Brownian path, 

r(r) = r + A 7 X(r), < r < q, (38) 

now parametrised by the (dimensionless) imaginary time r, < r < q. 
The path is specified by its position r in space, a particle species 7, 
a number of particles q, and a shape X(r) with X(0) = X(g) = 0. 
The positions of the q particles are located at points r(k — 1) on the 
path, k = 1, ...,q . The paths X r (r), r = 1, . . . , n, corresponding to n 
different loops are independent random variables 

(X?(t)X»(t')) x = 0, t + s (39) 

and identically distributed according to a normalized Gaussian mea- 
sure -D(X) with covariance 

(X?(t)X»(t')) x = J D(X)X»(t)X»(t') 

r ( T T> \ tt'i 

= 5^q min -,— ¥ , r = s, n, v = 1, 2, 3 . (40) 

L \q q J q z J 

The number q accounts for the quantum statistics of the species 7, it 
corresponds to grouping together q particles that are permuted accord- 
ing to a cyclic permutation of length q. The set of all possible loops 
(37) will be called the space of loops. It plays the role of an auxiliary 
classical-like phase space where methods of classical statistical me- 
chanics can be used. Note that for Bose or Fermi quantum statistics, 
the N particles are distributed into n loops C r , r = 1, . . . , n, according 
to their species and N = J2r=i1r- Maxwell-Boltzmann statistics are 
recovered if all (/-loops for q > 2 are disregarded. Then a loop C re- 
duces to a filament T and the covariance (40) reduces to (35) so that 
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in this case there is a one-to-one correspondence between filaments 
and particles. 

The generalisation of the Feynman-Kac-Ito formula to the many- 
body problem induces loop self-interactions and interactions between 
loops. The total energy of a system of n loops has three contributions: 

n 

U(£ r ) + U pot (£i,...,£n) + U A (£i,...,£ n ) . (41) 

r=l 

The potential energy J7 po t of n loops is the sum of pairwise interactions 
between loops plus the action of external potentials 

n n 

U^d, ...,C n ) = ^e 7r e 7s y c (£ r ,£ s ) + V cxt (£ r ) (42) 

r<s r=l 

where the interaction between two different loops is Coulombic 6 

VAAn = I* l ' " T ' - ^ |r(r)-rV)| - (43> 

Here, 5(f) = X^nL-oo ~~ n ) 1S ^ ne Dirac comb of period one, f = 
r mod 1. Hence V c (C r ,C s ) represents the sum of the interactions 
between the particles in the loop C r and the particles in the loop C s , 
and the factor 5(f—f') implements the quantum mechanical constraint 
of equal time interaction inherited from the Feynman-Kac-Ito formula. 
The term Ylr=i U O^r) ^ s ^ ne se ^ energy of the loops with 

U[C) - f f * f iT ' (1 - s ^ Wf - f, W)-Woi ■ (44) 

This is the sum of the mutual interactions of the particles within one 
loop. The factor (1 — £[ t ],[t'])> where [r] denotes the integer part of r, 
avoids counting the proper self-energies of the point particles; when 
q = 1, U(£) vanishes. Finally, 

U A (Ci,...,C n ) = -ifl fP-^i TdX r (r) • A(r r + A 7r X r (r),«(f)) 
r=l v P m lr c ■'O 

(45) 

6 A local regularization of the Coulomb potential has to be added when dealing with 
Maxwell-Boltzmann statistics. 
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where a(f) is the periodic extension of a(r), < r < 1 to all r. The 
phase factors in (45) arise from the interaction of the particles with 
the vector potential. They are the flux of the corresponding (periodic) 
magnetic field across the loops. 

The following remark is in order. In (43)-(45), r-integrals run 
from to q as a consequence of grouping together in a single path 
X(r), < r < q, all particles belonging to a permutation cycle of q 
elements (see [7], Chap. V, Section Al). Such integrals can as well be 
reduced to the interval < r < 1 by means of the identity 

/ dX(r)F(X(r),a(f)) = V / dX(r + m)F(X(r + m), a(r)) . 

J m=Q J ° 

(46) 

The notation in (43)-(45) is short and convenient. 

The total Gibbs weight on the space of loops (including the normal 
order constant Dn (20)) 



e- f3E>N exp 



-(5 U ( £ r) + U^d, ...,£„) + Uxid, £ n 

\r=l J _ 

(47) 

gives (up to normalisation) the joint probability distribution of n in- 
teracting loops in a realisation of the electromagnetic field having am- 
plitudes air). Individual loops have Gaussian weights defined by the 
covariance (40), thus calculations of averages on loops reduce in princi- 
ple to applications of the Wick theorem. One will also have to consider 
averages of stochastic integrals involving the line elements dX^ir). 
This is achieved by supplementing (39) and (40) by the expressions 

(dX?(T)X»(T>))x = ^(d^(r)X;(r'))x) dr, 



,8 l u,(e(T-T J )-'jjdT, for r + r' 

(48) 



WWW)x = S rs S, 



1 / d 



2 \dr 



1 



(^(r)X;(r)) x dr 



SrsS^u [ | - - )dr, for r = r' (49) 



15 



and 



(dX?(T)dX»(T'))x = (JLp ^(r)X;(r'))x) drdr' 

= (s(t - t') - -^j drdr' . (50) 

These formulae are in accordance with the middle point rule, which 
assigns the value 1/2 to 9{t — t')\ t=t > in (49) (see e.g. calculations in 
the appendix A of I). 

At this point we see that computations of thermal properties of the 
system of charges and field corresponding to the Hamiltonian (4) are 
entirely specified by the form of the Gibbs weight (47) on the space 
of loops together with the Gaussian distributions of the field ampli- 
tudes at(-) and loop shapes X r (-). Indeed, the Gibbs weight (47) is a 
functional of a(-) and X r (-), and Gaussian averages are uniquely char- 
acterized by the covariances (28), (29), (39), (40), (48), (49) and (50). 
Of course, calculation rules in the auxiliary space of loops have to 
be completed by appropriate formulae that relate quantities obtained 
in the loop formalism to the physical information of interest such as 
thermodynamic potentials or particle and field correlation. We shall 
not develop such formulae in general here but will present an appli- 
cation of this formalism to the determination of the particle density 
correlations in presence of the field in Section 6. 



5 The effective magnetic potential 

We are now in position to explicitly trace out the field degrees of 
freedom to obtain the representation of the matter statistical weight 
[e~^ HL ' R ] mSi t (31) on the space of loops. The corresponding distribu- 
tion is obtained by averaging (47) on the field variables, namely 

n 

e-f* D » exp[-/3 U{C r )\ exp[-/3C/ pot (£ 1 , . . . , C n )\ 

r=l 

x (exp[-PU A (£!,...,£„)] ) rad (51) 

From (45) and (32) one sees that exp[— f3U a(£i, ■ ■ ■ is a phase 

factor linear in the field amplitudes a^x^r) and aj^r). Since 

< • • • >rad is Gaussian, the average can be performed with the help of 
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the basic formula (written here for a single mode of the field) 
'exp \i f 1 dr(f(r)a*(r) + /*(r)a(r))l ) = 

\ l JO J / rad 

dr f dr'r(T)(a(T)a*(T% ad f(r') . 
Jo 



(52) 



To apply this formula we introduce the eigenmode expansion (32) of 
the vector potential in (45) 



-PU A (£i, ... ,£ n ) 



E (E jf dX r(r) ■ Kxir^j « kA (f) + c.c 



(53) 



where u£ a (t) collects the factors 



u kA (r) = P- 



-7r 



Application of the formula (52) gives 

(exp[-/ft7 A (£l, • • • , £n)])rad = 



(54) 



exp 



EE / dX r (r) ■ (ui x (r)y J2 

kA r=l ^ s=l 7 ° 



dX s (r')-K A (r'))C(A;,f-f') 
(55) 



We have used the fact that the covariance (27) is diagonal with respect 
to kA and C(k,r — t') is given by (28), (29). The remark made after 
(45) applies also here. In order to use (52), all r- integrals can as well 
be reduced to the interval < r < 1 by means of the formula (46). 
Then C(k,f — f') is the periodic continuation of C(k,T — r'), < 
t, t' < 1. Since u^ kA (r) = ±(u£. a (t))*, it is clear that by changing 
k — > — k, r — > s in (55) only the even part of C(k, f — r') contributes. 
One finds from (28) for r / 



Ceve„(fc,r) = -[C(k,r)+C(k,-r)} 

1 



n k cosh(/3^ k r) + 

cosh[/3na;k(|T| - 1/2)] 
sinh(/3^ k /2) 



(56) 
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whereas from (29) 



C cvcn (k,0) = n k , r = . (57) 
Introducing the explicit form of u kA (r) (54), equation (55) becomes 

(exp[-/3£/ A (£i, • • • , £n)])rad = 

exp (-p y / f^/^K ^-^cw 

x jf dX?{ T ) j*' dX;(r / )e ik '( A ^ x '-^- A ^ x ^ r '»C even (A;, f - f) ) • 

(58) 

The transverse delta function 5j^(k) results from the polarisation sum 

E^ = V-ns- = ^(k) • (59) 

A=l 

There is an important point to deal with before proceeding to 
the determination of the effective magnetic potential. The function 
C evcn (k,T) is continuous except for the point r = where it has the 
jump 

limC even (/c,r) - C cvcn {k,0) = - . (60) 

Although this point is of zero measure with respect to the Lebesgue 
measure, it cannot be disregarded when dealing with stochastic inte- 
grals. Indeed, when averaging over loops, the singular part S(t — r') 
in the covariance of stochastic differentials (50) will precisely select 
the value of C even (k,T — r') at r = t'. As an illustration, one can 
consider the X-average of (58) to linear order in the expansion of the 
exponential, namely 



dX?( T ) d X»( T ')e ik < X ^r(T)-\, s X s (T'))\ Ceven (k,f -f') . 

Jo Jo I X 

(61) 
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The average < • • • >x in (61) can be calculated by means of the Wick 
theorem, evaluating all contraction schemes. Contractions involving 
the product of stochastic differentials yield the term 

rf q '< dX^(r)dX^(r') >x / e ^(A 7r X r (r)-A 7s X s (r'))\ Ccvcn (k, f - ?) 

Jo Jo \ /x 

= S rs S^ r dr r dr' (5{t - r') - -) 

Jo Jo V Qr/ 

x / e ik.(A 7r X r (T)-A 7s X s (r'))\ Ceven (k, f - f) . (62) 



X 

In view of (60) the contribution of 5(t — t') in (62) is 

Srs^ u qrCcven(k,0) = 5 rs 5^ v lim C even (k,T - t') - 5 rs 5^ < ^- . (63) 

Then the contribution of the last term of (63) to the complete expres- 
sion (61) gives 

2irh 4. / i sr g 2 (k) 



n 
r=l 



R 3 S /, 
\ k 



(3j2Qrd lr =/3D N . (64) 



r=l 

The last line follows from the fact that we have n loops, each of them 
containing q r particles of species j r , so that Dn is the constant (20) 
arising from the normal order rule in the bosonic integral. At linear 
order, this constant exactly compensates the term —(3Dn occuring 
in the exponent of the total Gibbs weight (31). We conclude from 
this observation and from (62) that we can as well use the continuous 
extension of C cvcn (/c,r) (56) to r = and suppress the constant 
in (31), (51). A proof that this statement holds for all orders is given 
in Appendix A. 

We can now cast the field average (55) in final form 

/ Re 2 \ 
e-' 3D ^(exp[-/3C/ A (£ 1 , . . . , £ n )]) rad = J] exp I -^_>W m (£ r , C r ) 

(65) 

x exp -PYl 

W m (£ r ,£ s ) 

V r<s J 

Here we have introduced the effective magnetic potential 

x rd^(T)e- ik - A ^ x '- (T) r'dXZ^e^'^W C(k)Q(fe, f 

Jo Jo k ' 
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To obtain (65) and (66), we have separated in (58) the terms r = s 
refering to the self energies of loops from the terms r ^ s giving rise 
to pairwise loop interactions. The function 

a ' fc ' T » = 2s m hX^) C0Sh[Vt(|r| - 1/2)1 
A ph fe\ eW'(M-i) + e -A ph fc|r| 

—) ' M - 1 (67) 

is, up to the factor A p h&, the even part (56) of the covariance of the 
free photon field written in terms of the photon thermal wave length 
A p h = Phc. In view of the discussion following (60) and the result 
of Appendix A, it is understood that this function is given by the 
formula (67) including the point r = and the factor e~@ DN has 
been cancelled in the right hand side of (65). The r-periodic function 
Q(k,f), Q(k, 0) = Q(k, 1), is normalized in such a way that it equals 
one when the electromagnetic field is classical 

lim Q(k,r) = 1 . (68) 

A ph ^0 

In this limit, the magnetic potential W m (£ r ,£ s ) reduces to formula 
(82) of I where radiation has been treated classically. Hence all effects 
due to the quantum nature of the photon field are contained in the 
sole function Q(k,r). 

The Gaussian integration of the radiation field has provided the 
sum of pair potentials (65) between loops as in Paper I. Then, ther- 
mal averages of particle observables calculated with the normalized 
reduced density matrix pl,r (10) have a simple structure when ex- 
pressed in the system of loops. Combining (51) and (65), one forms 
the complete effective Gibbs weight (up to normalisation) 



exp 



x exp 



\r=l ~ / 



n 

- /3(V P ot(£i, ...,£„) + £ e 7r e 7s >V m (£ r , £,))] (69) 



r<s 



comprising one-loop and two-loop interactions. This structure allows 
the use of standard diagrammatic methods of classical statistical me- 
chanics, like Mayer graph expansions. This is illustrated in the next 
section, where large-distance asymptotic particle correlations are in- 
vestigated. 
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Note that as in Paper I, it is unlikely that pl,r can be cast in a con- 
venient operator form pi,R oc e~^ HeS ^ Pi ' ri ^ depending on the original 
quantum-mechanical momenta and positions {pi,i"j} of the particles. 
Again, the magnetic interaction W m (66) is a two-times functional of 
the Brownian loops reflecting the photonic bath environment. It lacks 
the equal-time constraint necessary to come back to a simple operator 
form by using the Feynman-Kac-Ito formula backwards [15]. 



6 Asymptotic particle correlations 

We determine the behaviour of the particle density correlation in 
presence of the thermalized quantum electromagnetic field in the two 
regimes (2) and (3) discussed in the introduction. 



6.1 Partial screening of the Coulomb interaction 
by thermal photons in the range Amat^Aph^r 

In the regime (2), r is larger than any typical length of the model. The 
asymptotic analysis of the correlation is based on the large-distance 
behaviour of the part of the interaction 7 

W{T a ,T b ) = W c (Fa,Fb)+Wm(Fa,Fb), (70) 

which is responsible for the power-law decay. In this formula, W c (J- a , J-},) 
is the residual interaction (due to quantum fluctuations) that is left 
when Coulomb divergencies are resummed in Mayer graphs (see for- 
mula (28) of I). It has the asymptotic dipolar form 

W c {F a ,Fb) ~ |r - r b | -> oo 

1 



Jo Jo 



ds b (5(s a -s b )-l) (\ la £ a (s a ) ■ V Ta ) (\ lb £ b (s b ) ■ V rb ) 

I 1 a r fe | 

(71) 

It turns out that the large-distance asymptotics of W m (/ a ,fj), de- 
termined by the small- /s behaviour of the integrand of (66), are also 
dipolar. Indeed, we first observe that Q(k, r) is an analytic function 
of k and has the small- expansion 



B(lt , T ) = i + 2s«! 



1 

- |r| + ~ 6 



+ 0((X ph k) 4 ) . (72) 



7 Exchange effects are short ranged and play no role here. Only one particle loops, i.e. 
filaments, are considered. 
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Inserting this in (66) gives 

Wm(Fa,Fb)~W m (f a ,T b ) 

x /"dtf (r) fde b (r f ) 
Jo JO 



j3yjm a m b c 2 



r dk 



3 ik-(r a -r 6 ) 



(r - r 



l\2 



T ~ T + 



fc 2 

|r a - r b | 



— ► oo 
(73) 



The first term in the r.h.s of (72) leads back to the effective magnetic 
potential W m associated to the classical electromagnetic field (formula 
(22) of I). In the second term, the k~ 2 factor in the integrand of (66) 
has been cancelled by the term of second order in k of (72) and we 
have set k = in the exponentials of the paths £ a ( T ) an d £fc( r ')- in 
this way, we have retained the lowest order singular part in k in the 
last term of (73). This part is — k^k u /k 2 coming from the transverse 
delta function (59). The double stochastic integral in (73) is calculated 
with the result 



/dr/dr' ^(r)^(r') 
Jo Jo 



d 2 



drdr' 



(r - r') 



/\2 



I „ 1 

r-r ' + - 
6 



2 /dr/dr' (S(r-r')-l) 
Jo Jo 



'0 Jo 
In virtue of the identity 



A P h 



(3^/m a m b c 2 



A a A& 



equation (73) eventually reads 



-A a A 



dr/dr' («5(r-r')-l) ^(r)^(r') 
o Jo 



dk 



ik-(r a -r b ) 



(74) 



(75) 



(76) 
4vr£^ 



Performing the Fourier transform, we see that, up to the sign, the 
second term in the r.h.s. of (76) is identical to the asymptotic tail 
(71) of W c . The latter is therefore exactly cancelled in the total in- 
teraction YJ{T a ,T b ) = WciFaifb) + ^(fs/t) as |r a - r 6 | -> oo. 
We conclude that in the region r S> A p h the dominant part of this 
algebraic Coulombic tail is screened by thermalized photons. The tail 
of the interaction 



(77) 
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reduces therefore to the pure unscreened effective magnetic current- 
current interaction W m (J- a ,J-b) induced by the classical field, whose 
asymptotic dipolar form is given by formula (25) of I. 

We can now follow the asymptotic analysis presented in Section 
V of Paper I to show that the tail of the correlation exhibits again a 
generic r~ 6 decay. All statements made there regarding the magnetic 
potential with the classical field W m hold for the magnetic potential 
with the quantum field W m . The transversality argument used to 
show the vanishing of the convolution element (1.50) works identically 
provided that the rotationally invariant function Q(k, S2 — fife) (67) is 
included in the definition of the tensor T U2 (k, si, s&) (1.51). This tensor 
still transforms in a covariant manner under rotations of k, so that its 
contraction with the transverse delta function cancels (1.50). Similar 
modifications done in the other convolution elements mentionned after 
Eq. (1.51) imply that W m does not contribute to the W-convolution 
chains occurring in (1.49). The dipolar character of the large-distance 
interaction W then ensures that the correlation function decays as r~ 6 . 
However, the amplitude of this decay is now affected by the partial 
screening (77) which is due to the quantum nature of the photonic 
bath. 

In order to illustrate this point, let us determine the coefficient of 
the r -6 decay at lowest order in h. Proceeding word for word as in 
Section V of Paper I, one sees that this decay is eventually governed 
by 

i[-/?e 7l e 72 W(^i,^ 2 )] 2 (78) 

with root points dressed by classical correlations and evaluated at 
lowest order in h. Since W = W c + W m depends on h solely through 
the couplings A mat /c in W c , and A ma t/c, X p hk in W m , evaluating these 
potentials at lowest order in h amounts exactly to selecting their large- 
distance {k — > 0) asymptotic behaviour. The Coulombic dipolar tail 
of W c is therefore cancelled by the photon-induced partial screening 
(77), and the large distance behaviour of the two-particle truncated 
correlation in the semi-classical regime (high-temperature or lowest 
order in h) reads: 

PT(7a,r a ,7 b ,r b ) — ^ dr n i(la, 7i, r) /dr nfe, 7 6 , r) 

71,72 \- J - 1 \- J 

e 2 e 2 1 

x r V 2 r J -i6- ( 79 ) 
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This corresponds to omit the Coulombic part of the correlation calcu- 
lated in I, Formula (53). Only the current-current interaction induced 
by the thermal motion of the particles contributes to the tail (79) in 
the regime r S> A p h- 

The analysis of the particle-charge and charge-charge correlation 
function can be performed in the same way. As recalled in (1.46) and 
(1.47), the Mayer bonds are built from a rapidly decaying resummed 
potential $ e iec and the quantum asymptotically dipolar potential W 
(70). When the charge observable is considered, following the dressing 
method described in Sect. VI. A. 3 of [7], one sees that an additional 
screening factor involving <I> c i e c occurs at the root points of the graphs. 
This generically weakens the decay of the particle-charge correlation 
to r~ 8 and that of the charge-charge correlation to r~ 10 . As for the 
particle-particle correlation (79), the amplitudes of the tails are again 
determined by W 2 . Because of the asymptotic cancellation of the 
Coulombic part in W (see (77)), these amplitudes also inherit the 
small relativistic factor (f3mc 2 )~ 2 . 



6.2 Predominance of electrostatic correlations 
in the range A mat «r« A p h 

Let us now focus on the second regime, A mat <C r <C A p h, i.e. we 
consider the correlation between particles that are separated by dis- 
tances much smaller than the wavelength of thermalized photons. We 
first give a rough estimate of the order of magnitude of W m (f s ,f(,) 
relative to W c (J r a , ^5). In this aim it is convenient to scale the Fourier 
variable as k — ► k/r, r = r a — r b , yielding in (66) 

1 1 f dk 

I3yfm a m b c 2 r J (2ir) 

(80) 

Since A p h/r is now a large number, it is not allowed to expand 
Q(k/r,r — t') for small k, but we note from (67) that this function 
is of the form A p h/r times a bounded function of A p h/r. Therefore 
Q(k/r, t — t') cannot grow faster than A p h/r. Then the order of mag- 
nitude of W m is at most 

W m = , — i= 2 O ft*) . (81) 
/jy/m a m b c z r \ r J 



W m (f a ,f b ) = ^=3: /7^e ik ^de(r)e ik ^/ r ^W 



r = - 

r 
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On the other hand, one sees from (71) that the order of magnitude of 
W c for r > A mat is 



(82) 



Combining (81) and (82) together with (75) gives 



.Xph 

Hence, in the range (3), the total interaction 



W m = W c & | ~~ J (83) 



W = W c + W m = W c ( 1 + O [ ^- ) ) (84) 



ph 



is given by its Coulombic part up to a small correction. It is there- 
fore expected that all predictions on correlation decays are the same 
as those derived from pure electrostatics up to terms that vanish as 
r/Xph —> 0. This reasoning is mathematically not complete since when 
(80) is used as a bond in Mayer graphs, loop averages and wave num- 
ber Fourier integrals have to be performed first and shown to yield 
finite values. As an example we establish in Appendix B the precise 
estimate 

«(^)>^ ~ 120A(W^ a ,F b ))^ b (^) 3 , A < oo 

(85) 

as A cut /r -> 0, X a /r, X b /r -> 0, A ph /r -> oo, implying 

{wl{^T b )) iM = {wl{F a ,F b )) iM ((^) 3 ) (86) 

in the range (3). Thus the square fluctuation of W m (entering e.g. in 
(78) for the evaluation of the correlation) is negligible compared to 
that of W c - 



7 Concluding remarks 

In this paper, we have presented a formalism adapted to the study of 
non relativistic matter in thermal equilibrium with the photon field. 
In the joint functional integral representation of matter and field, the 
field variables can be integrated out, yielding an effective classical- like 
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statistical description of the state of matter. As a first application, 
we have shown that the cloud of thermalized photons participates 
in the screening of the Coulomb potential by supressing the dipolar 
electric contribution to the r -6 tail of the particle correlations, as il- 
lustrated in (79). In electrolytes at room temperature, both the de 
Broglie and classical Debye lengths are in the range of 10 _10 m (a few 
Angstroms). Moreover, the parameter \J f3mc 2 is of order ~ 10 5 so 
that A p h ~ lCT 5 m, see (1). It is known from [20] that the crossover be- 
tween Debye-Hiickel (exponential) screening and quantum (algebraic) 
screening occurs at distances of about 60 times the Debye screen- 
ing length. Consequently, the further reduction of the amplitude of 
the correlation tail by photon screening occurs at even much larger 
distances and with an exceedingly small amplitude. This makes the 
phenomenon probably hardly observable in such systems. At the con- 
ceptual level, it is however an interesting effect of the thermal radiation 
that, to our knowledge, has not been exhibited in the literature before. 

The effective magnetic potential W m (66) defined in Section 5 em- 
bodies in an exact manner orbital diamagnetic interactions, namely 
interactions between currents due to thermal motion of charges. This 
current-curent interaction is at the origin of the correlation tail (79). 
The order of magnitude of W m is by a factor (/Jmc 2 )" 1 smaller than 
that of the electrostatic potential. One should however be aware that 
W m is not the unique source of relativistic effects. A preliminary in- 
vestigation [21] shows that the Pauli coupling terms of spins with the 
field contribute to the correlation tail at the same order (j3mc 2 )~ 2 as 
that found in (79) as a consequence of pure orbital magnetism. More- 
over, the non relativistic form of the particles kinetic energy in the 
Hamiltonian (4) has itself c~ 2 corrections (e.g., spin-orbit interaction, 
Darwin term) that will likely contribute to the asymptotic form of the 
particle correlations. Hence a complete determination of the particle 
correlations tails at order (/3mc 2 )" 2 will require further investigations. 

The tools developed in this paper lend themselves to a detailed 
microscopic study of various problems. In view of the brief report [22] 
questioning the findings of paper I on electromagnetic fluctuations, we 
aim to revisit the problem in the case of a quantized field. Thermal 
broadening of spectral lines and retardation effects on van der Waals 
forces between recombined atoms or molecules in a medium at finite 
density and temperature could conveniently be studied within this 
formalism. Indeed, the latter situations involve quantum mechan- 
ical binding which is not perturbative in the matter-field coupling 
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constant. Standard many-body Feynman diagram techniques would 
necessitate infinite resummations to describe bound state formation, 
whereas cluster expansions in the form presented in [23] (properly 
generalized to include the full electromagnetic coupling) give a direct 
access to recombined entities together with their interaction with the 
radiation field. Finally the theory of the Casimir effect has received 
much attention recently. It is now conceivable to elaborate a full mi- 
croscopic theory of this effect by extending the analysis presented in 
[24] to TQED. We plan to address these questions in future works. 

Appendix A 

As seen in the first order calculation leading to (64), the compen- 
sation of the constant Dn comes from the particle self energies r = s. 
It is therefore appropriate to single out in the exponent of (58) a di- 
agonal r = s term and write its X r average as (dropping now the 
particle index r) 



/ = ^exp 



r j 3 w i' q r q 



F(X)\ , (87) 

/ X 



x e -ik.A(X(T)-X(T' )lVti/(fc)f _ fl) 

where we have set for brevity 

T, u (k, r) = ' 9 l ik) ^(k)C even (fc, r) . (88) 

In (87) F(X) is a functional of X containing all possible other depen- 
dences of X in (58). Expanding the exponential in (87) gives 



n ! \ 1J -y (2vr) 

71=1 \j = l J K ' 



f[ dX* (Tj) J* dX v > (^) e -*j-*(Xfa)-X(^) j F(X) 



X 



II I/-. ] . (89) 



We call a matched contraction the contraction of a pair of stochastic 
differentials < dX* (Tj)dX^ (rj) >x= S^ jVj (Sfa - rj) - d^dr', 
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where times have the same index j. It is clear that the 5{rj — rj) 
occuring in matched contractions will evaluate r HjVj{kj,fj — fj) at 
Tj = Tj. Such matched contractions can only arise from the product 
Yl] =1 dXM (Tj)dX u i (rj) in (89). Contraction between a differential 
from this product with a differential occuring in -F(X), or contrac- 
tions within -F(X), will always involve two time arguments belonging 
to different C eV en functions. They are of the type 

(J* dX»( T ) J* dX v {a)T^ (k, f - f')T uv ,{k\ a - a')^ = 

J* dT J* da (8{t - a) - -^j T^,(k, f - f ')r> (k', a-a') = 

q [ dTT^ik^-f^T^ik' ',t-ct')- f drT w '(A;,r) f daT ^ (k' , a) 
Jo Jo Jo 

(90) 

For such contractions, C eV en(^ ) 7") can be treated as a continuous func- 
tion everywhere since in integrals of the type (90) the discontinuity 
(60) at the single point r = is irrelevant. 

To evaluate the X average in (89), we select therefore terms having 
exactly m matched contractions, < m < n. Because of the invari- 
ance of the product under exchange of its factors there are 
n\jm\{n — m) ! such terms giving the same contribution. This leads 
to 



n ! 



n=l 
m 



m=0 



m !(n — m) ! 



dkj 




So ^ So dT ' jF ^ kj,fj ~ fj '^ \ ^ Tj ~ ~ \) 

(B(X)) n - m F(X)\ 



e -ik,-A(X(r,)-X(rj)) 



unmatched 

(91) 



with 



B(X) = 
dk 



SwfSo dx " (T) f dx "(r>- ik ' mT) ~ MT ' ))T ^ f ~ f ') 



(92) 
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The square bracket in (91) is the result of m matched contractions. 
In the average < • • • > unmatched) ah matched contractions are omitted. 
In (91), we further expand the product of (j)(Tj — rj) — 1/q^j factors 
and perform the 5 function integrations leading to 



oo 


-13)" 
n ! 


n=l 


"{"J 


r dk 


(2tt) : 



II I lit I 

En ! v - m ! 

m \(m — n) ! I \(m — I) ! 

m=0 v J 1=0 y ' 

jF^k, 0)) ' ((B(X)) n - m (D(X)r~ l F(X) 



unmatched 

(93) 



with 



£(X)=-- ll^Fdr [ 9 dr'e-^W-Mr')) r^.f-f ) . 
qj (Ztt) 6 Jo J 

(94) 



Finally, rearranging the sums yields 
dk 



I = exp -pq 



IW(M)) (e- B ^e- D ^F(X) 



[ZTT)^ J \ / unmatched 

(95) 



It is seen from the definition (88) that f3q J f0yzF^(k,O) is equal 
to the constant — qd plus the contribution of C even (k,T) extended by 
continuity at r = 0, exactly as in (62)-(64). Since the loop shapes 
X r , r = 1, . . . , n, are independent random variables, the same calcu- 
lation can successively be carried out for n loops, providing a factor 
e /3D N that cancels the factor e _/3Djv due to normal ordering in (51). 
Performing the procedure (87)- (95) backwards after this cancellation 
thus shows the validity of formula (65), where the effective magnetic 
potential W m (66) is defined with the continuous function Q(k, r) (67) 
for all r. 

Appendix B 

The £ a ,£ft average of W m reads in terms of the scaled variables 
q a = k a r and q b = k b r 

/V0 2 (T T\\ = - Lpf^i ^1 ^ ^Eii 

(96) 
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where we have introduced the function of dimensionless parameters 



F ( Acut K h Aph\ = f d Qi f d( l2 
V r ' r ' r ' r J J qi< ^_ (2vr)3 J (2vr)3 

^¥c(Qi)^(g 2 ) (j«<¥e*<T)-$s.W) e uh-(¥<.w-^(0) 



o JO JO Jo / a *. v r / \ r / 



(97) 

The ultra-violet cut-off functions and gfe) have been replaced 

by the appropriate restrictions of the domains of integration. Then 

Aph ^ _ j. ^ / A cut A a Aft Aph 



/ A cut a„ \ 1 r ' r ' , 



J (2vr)3 7 (2tt)3 q \q\ VWi^WW X 

Vo Jo Jo Jo / £ o ,£ 6 v r / V r / 

(98) 

The £ a , £ fe average is evaluated according to the rule (50). Using the r 
periodicity of the function Q(k,r), the square bracket in (98) becomes 

(A ph /V)V<72 



2(1 — e -( A phA)9l)(l _ e -(A p h/r)92) 

(1 _ g-(A p h/V)(«l+<22) g-(A p h/r)gi _ e -(A ph /r)(?2 \ 
(A P hA)(9i + 52) (A P hA)(gi - <?2) y 

x , A p h gi<72 Aph , , 

~6vx8 v s — — ( — — -rr, >oo 99) 

as shown by an explicit calculation of the r integrals. Inserting (99) 
in (98) and performing the vector sums leads to 

^) ~ ^ (100) 
r I r 
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with 



fdq 1 f 



A= /^- I dq2 e^+rt* ^ ( (9i ; 9 2)2 - 3] - 

J 2( qi +q 2 

(101 



(2tt) 3 qiq 2 \ q\q\ ) 2(gi + q 2 



< 00. 



Introducing the representation 1/(51 + q 2 ) = dte~ t ^ qi+q2 \ the qi 
and q 2 integrals can be performed independently and each of them 
behaves as t~~ 2 as t — ► 00, assuring the convergence of the i-integral. 
Since (W^a, F b )\ ~ A 2 Af/120r 6 as A a /r, \ b /r -> 0, one obtains 
(85). 
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